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. Abstract 

o . 

' Let T^^^ be the first hitting time of the point 1 by the Bessel process with index G M 

^ ' starting from x > 1. Using an integral formula for the density qi^\t) of T^^\ obtained in 

' [5], we prove sharp estimates of the density of T^^^ which exhibit the dependence both on 

^sO ■ time and space variables. Our result provides optimal uniform estimates for the density 

^ of the hitting time of the unit ball by the Brownian motion in M", which improve existing 

I bounds. Another application is to provide sharp estimates for the Poisson kernel for 

Ph ' half-spaces for hyperbolic Brownian motion in real hyperbolic spaces. 

g ' 1 Introduction 

I Bessel processes play a prominent role both in the theory of Brownian motion (see [15] and 
I> ' [H]) as well as in various theoretical and practical applications. The n-dimensional Bessel 
^ I process appears quite naturally as the n-dimensional Euclidean norm of Brownian motion; the 
i/^ ', more intriguing applications are related to the celebrated Ray-Knight theorems describing the 
^ i behaviour of the local time of Brownian motion in terms of two-dimensional (quadratic) Bessel 
0^ ' process (see [18], [20] or [2]). There is also an intimate relation between Bessel processes and the 
^ ' geometric Brownian motion [12] (the so-called Lamperti's theorem - see Preliminaries). Bessel 
^ ■ processes also appear when representing important jump Levy processes by means of traces of 
>- . some multidimensional diffusions [Bessel-Brownian diffusions); see [IH], [1]. Another important 
^ i application consists of the fact that the hyperbolic Brownian motion, i.e. the canonical diffusion 
^ [ in the real hyperbolic space, can be represented as subordination of the standard Brownian 
■ - - ' motion via an exponential functional of geometric Brownian motion. Thus, from Lamperti's 
theorem, the Poisson kernel of half-spaces for hyperbolic Brownian motion can be represented 
via subordination of the standard Brownian motion by the hitting time of a Bessel process 
(see [5]). We exploit this relationship for obtaining the precise bounds of the Poisson kernel of 
half-spaces for hyperbolic Brownian motion with drift. 

Our main goal is to study estimates of the density of the distribution of the first hitting 
time Ta'^^ of a level a > by a Bessel process starting from x > a. Our approach is based on 
an integral formula for the density of Ta^\ given in the paper [3] (see also [3]). This formula. 
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although quite complex, proves to be a very effective one. We deal with processes with non- 
positive indices, however our conclusions are valid also for positive ones, since the processes are 
equivalent on the sigma- algebra generated by the process before hitting 0. 

To formulate our results denote by R^'^^ = {R['^\t ^ 0} a Bessel process with index /i e M. 
We note by P^f^-* the probability law of a Bessel process i?*^^^ with index n on the canonical 
paths space with starting point Rq^^ = x, where x > 0. Let us denote the first hitting time of 



•-0 

the level a > by a Bessel process with index fi 



Ti'^)=inf{t>0;i?r^ = a}. 

Our main result is a sharp estimate of the density of the hitting distribution if x > a. By the 
scaling property of Bessel processes it is enough to consider x > a = 1. The distribution of t/'*'' 
is well known for /i = ±1/2 when both distributions have the 1/2-stable positive distribution 
with scale parameter A = x — 1, however incomplete for /i = 1/2. Our upper and lower estimates 
are comparable in both time and space domains. We have the following: 

Uniform estimate for a density function of T^^K 

For every x > 1 and t > we have 

Pi^\Tl''^ e dt) _ f 1 \ e-("-i)V2t -^2|M|-i 

It -{x-1) y^^^2t.J tH-i/2 + a;lMl-i/2' 

Moreover, we have 

Pi°)(Tr G dt) ^ ,,^_(._i)V2t(^ + ^ 1 + 



1 e 



dt ' ' xt^l'' (l + log(l + |))(l + log(t + x))' 

Here f ^ g means that there exist strictly positive constants ci and C2 depending only on 
such that ci ^ f /g ^ C2. 

Our setup includes the case of the hitting distribution of a unit ball by a Brownian motion 
starting from the exterior of the ball. To state this result, let a'-"'-* be the first hitting time of a 
unit ball by n-dimensional Brownian motion = {W^'^\t ^ 0}, i.e. 

a(") = inf{t>0;|iy/"^| = l}. 
Uniform estimate for hitting the unit ball by Brownian motion. 

For Wq^"^ = X G M" such that |x| > 1 we have 

/2t -| 

n > 2, 



P^'(f^^"^ e dt) _ 


X 




- 1 e-{l^'l-i)V2i I 






dt 


X 


f3/2 t(n-3)/2 _^ 


X 


{n-3)/2 



for every t > 0. Moreover, we have 

-{\x\-lf/2t 



PHa^^^ G dt) 


X 




- 1 


dt 


X 







+ 1 + log 


X 






t3/2 (l+log(l+ *))(1 


+ 


log(t + 


\A)) 



To the best of our knowledge even for the planar Brownian motion our results are new and 
considerably complement existing results obtained in [9] (see also [6]), where the estimates are 
only sharp in the region t ^ |xp with sufficiently large starting point x. If t < |xp the bound 
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obtained in |9] has an exponential term of the form exp{—Ci-^} with different constants Ci, C2 
for the lower and the upper estimate, respectively. We remove this obstacle and provide sharp 
estimates which are of the same order in the full range of t and x. We also provide sharp 
estimates for the survival probability Pl^\t < t[^^ < 00). The asymptotic result if t — cxd 
is due to Hunt [10] in the case of the planar Brownian motion and Port [17] in the context of 
Brownian motion in higher dimensions. Recently, a result about the asymtotic behaviour of 
the hitting density for the planar Brownian, when t 00, was established by Uchiyama [T^ . 
His result gives a very accurate expansion of the hitting density provided |xp is small relative 
to t, that is when the exponential term is negligable. When this is not true the error term of 
the expansion in [19] may be much bigger than the leading term. Therefore our estimate in the 
case of the planar Brownian motion is much more accurate, when we assume that t is not too 
large with respect to since the impact of of the exponential term may be significant and 
it is reflected in our estimates. Moreover, we give a very exact estimate of the density in the 
situation when t is small relative to |x| (see Lemma [Hand Remark [1]). 

The organization of the paper is as follows. After Preliminaries, in Section 3, we provide 
uniform estimates of the density function of the first hitting time T^^\ This section is basic 
for further applications, which are collected in the next section. We first provide the estimates 
of the survival times of a killed Bessel process and, finally, compute the precise bounds of the 
Poisson kernel of a halfspace for hyperbolic Brownian motion with drift. Appendix contains 
various estimates of quantities involved in the basic formula for the density function of T^'^\ 
which, although quite laborious, but at the same time, are indispensable ingredients of the 
proof of the main result. 

Throughout the whole paper f ^ g means that there exists a strictly positive constant c 
depending only on fi such that c~^g ^ f ^ eg. If the comparability constant will also depend 
on some other parameters 71, 72, . . . we will write f g, e = 0(71, 72, . . . ). Also in a string of 
inequalities a constant may change from line to line which might not be reflected in notation. 
Moreover we consistently do not exhibit dependence of c on in inequalities of type f ^ eg. 



2 Preliminaries 



2.1 Modified Bessel functions 

Various formulas appearing throughout the paper are expressed in terms of modified Bessel 
funetions and K^. For convenience we collect here basic information about these functions. 

The modified Bessel functions of the first and second kind are independent solutions to the 
modified Bessel differential equation 

z^y" + zy' - {^^ + z^)y = 0, 

where 'i? G M. The Wronskian of the pair {K^{z), I^{z)} is equal to 

W{K^{z),h{z)] = l4z)K^+i{z) + h+i{z)K4z) = ^. (1) 

In the sequel we will use the asymptotic behavior of and at zero as well as at infinity. 
For every -i? ^ we have (see [1] 9.6.7 and 9.6.12) 

I4r) = + 0(r-^+^) , r^0+,^>0. (2) 
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For ?9 > 0, we have ([I] 9.6.9 and 9.6.13) 

Mr) = r^0\ (3) 

where g{r) = f{r) means that the ratio of g and / tends to 1. Moreover, in the case "(9 = 0, we 
have (see P 9.6.13) 

J^o(r) = -log^Jo(r) + 0(l), r ^ 0+. (4) 
The behavior of and at infinity is described as follows (see [1] 9.7.1, 9.7.2) 

/,(r) = ^(l + 0(l/r)), r^oo, (5) 

K4z) = ^e-\l + 0{l/z)), 1^1 ->oo, (6) 
where the last equality is true for every complex z such that |ar(72;| < |7r. 

2.2 Bessel process and exponential functionals of Brownian motion 

In the following section we introduce notation and basic facts about Bessel processes. We follow 
the exposition given in [13] and [H], where we refer the Reader for more details and deeper 
insight into the subject (see also [18]). 

We denote by P^^'' the probability law of a Bessel process R^'^^ with index fi on the canonical 
paths space with starting point R^^^ = x, where x > 0. Let J^^^^ = a{Ri^\s ^ t} be the 
filtration of the coordinate process r[^\ The state space of R^^^ depends on the value of /i and 
the boundary condition at zero. For simplicity, in the case — 1 < yU < (then the point is 
non-singular), we impose killing condition on 0. However, the exact boundary condition at 
is irrelevant from our point of view, because we will only consider the process i?^^-* up to the 
first hitting time of the strictly positive level. 

Let us denote the first hitting time of the level a > by a Bessel process with index fi 

Ti^) = M{t>0;Rl^^ =a}. 

Observe that for /i ^ we have Ta^^ < oo a.s. and P^^^Ta^^ = oo) > whenever /i > 0. Using 
the scaling property of Bessel processes, which is exactly the same as the scaling property of 
one- dimensional Brownian motion, we get for every b > and t > 

Pt\Tl;^^ <t) = Pi'\b'Ti^^ <t), X > a > 0. 

Therefore, from now on we do assume that a = 1 and x > 1. We denote the density function 
of T^^^ with respect to Lebesgue measure by qi^\ i.e. 

gi^^W = ^^^^^^^1^, t>0,a;>L 

We have the absolute continuity property for the laws of the Bessel processes with different 
indices 
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Here Tq^^ denotes the first hitting time of by R^^^K If z/ ^ then the condition {Tq'^^ > t} can 
be omitted. In particular, for z/ = — /i, where /x ^ we have 



r(-M) 



(^)'', Pi"-) - a.s. on {t(-'') > t}. 



<li'\t)=(^y%i-^\t). (8) 



Consequently, for /i ^ and a; > 1 we get that 

X J 

We denote by i? = {Bt,t ^ 0} the one- dimensional Brownian motion starting from and 
by B^^^ = = Bt + fJ,t,t ^ 0} the Brownian motion with constant drift /i G M. The process 

— {x exp(i?j^''), t ^ 0} is called a geometric Brownian motion or exponential Brownian 
motion with drift /i G M starting from a; > 0. 

For X > 1 let r be the first exit time of the geometric Brownian motion with drift /i from 
the set (1, oo) 

r = inf{s > 0; X exp(i?s + fis) = 1}. 
We have r < oo a.s. whenever /i ^ since then inft^o^^^H^) = 

— oo. 

For X > we consider the integral functional 



Ai^\t)= [ = /" exp{2B, + 2fis)ds. 

Jo Jo 



The crucial fact which establishes the relation between Bessel processes, the integral func- 
nal Ai^^ and the geometric Bi 
a Bessel process R^^^ such that 



tional Ai^^ and the geometric Brownian motion is the Lamperti relation saying that there exists 



Consequently, we get 



xexp(i?l^^)=i?^^^^^, t^O. 



Ai^\T)^T['\ (9) 



2.3 Representation of hitting time density function 

We recall the result of [5], where the integral formula for the density of A^x was given. 

According to ([9]), as immediate consequence, we obtain the formula for ^\t)- Note also 
that in the paper [5] different normalizations of Brownian motion and different definition of 
geometric Brownian motion were used and consequently we have gi ^\t) = g^(t/2)/2, where 
qfj.(t) is the density function considered in [5]. 

Theorem 1 ([Byczkowski, Ryznar 2006]). For fi ^ there is a function Wx such that 

<li-'\t) = X'-^ + [e-^/^^ - 1) w,{v)dvj , (10) 

where n = k{v) = (A + vY — = v{2X + v), and X = x — 1. 
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The function w\ appearing in the formulas is described in terms of the modified Bessel 
functions and J^. The function K^{z) extends to an entire function when /i — 1/2 is 
an integer and has a holomorphic extension to C \ (— oo,0] when /i — 1/2 is not an integer. 
Denote the set of zeros of the function K^{z) hj Z = {zi, Zk^^} (cf. [7], p. 62). Recall that 

= H — 1/2 when — 1/2 G N. For /i — 1/2 ^ N, /c^ is the even number closest to /i — 1/2. 
The functions and K^_i have no common zeros. 

The function w\ is defined as a sum of two functions 



Wx{v) = Wi^x{v) + W2,x{v), (11) 



where 



and 



Moreover the moments of n with respect to w\{v)dv can be computed in the following way 

x'^-i/^^/i^ - l/4)/2x = / wx{v)dv, (12) 

and, for /i > 1/2, we have 

2x^-^/2 = / Kwx{v)dv. (13) 







We use also the following representation of qx (t) for ^ 1/2 (see [5] (24)) 

-A2/2t 



where / = [/x + 1/2] if /i ^ N, and 1 = ^—1/2 otherwise. 

3 Uniform estimates of hitting time density function 

Throughout of the rest of the paper we denote A = x — 1. The main result of the paper is the 
following uniform estimate for a density function of T[^\ 

Theorem 2. For every x > 1 and t > we have 



. I \ g-AV2i ^2|m|-1 

^i^Ht)^A(^^J ,3/. ,H-i/2 + ,H-i/2 > /^^O- (15) 



Moreover, we have 



^ g A /2i 1 + log X ^ 2 

9™ (i) <! ^J^Jl + 'ogf)(l + log()' ' > 

-^3,1/2^3/2' t ^ 2x 



or equivalently 

^1/2 l + logx 



(0) 



it) ^ Ae-V- Oi+I 



xtV^ (l + log(i + i))(i + log(t + x))' 



As corollary, putting fi = n/2 — 1, we get the corresponding result for n- dimensional Brow- 
nian motion. 

Theorem 3. Let a^"'^ be the first hitting time of a unit ball by n-dimensional Brownian motion 
l^W = ^ 0}, I.e. 

cr(") = inf{t > 0; = 1}. 

Then, for ' = x such that \x\ > \ we have 



P^(ct(") e dt) _ 


\x\ 




. I ^-{\x\-lf /2t I 






dt 


\x\ 


^3/2 t(n-3)/2 _^ 


\x\ 


(n-3)/2 



n>2 



for every t > 0. Moreover, we have 



PHa^"^^ e dt) 


\x\ 






\+t)^/^ 1 + log 


\x\ 




dt 


\x\ 




t3/2 (l + log(l + ^))(l+log(t + 





The proof of the main theorem follows from Lemmas El [71 [H [9] given below. We use 
the crucial estimates of the function wx{v) and its components wi^x{v) and W2,\{v) given in 
Appendix (see Lemmas [HI [16] and [TTj). 

The following lemma provides satisfactory estimates of q^x^^ (t) in the case when t is small 
relative to x. 

Lemma 4. We have the following expansion 

where the error term satisfies the following estimate 

\E{t,x)\^C-{ViAU. 

X A 

Moreover, for ^ fi < 1/2 we have 
for every x > l,t > 0. 

Proof. By the basic formula ([TOj) . with application of (fT2|) . we obtain 



xi"- 




'0 


1 - 




t 




8 


X 


1 - 


4/i2 


t 
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Observe that, by Lemmas fT4l [T6] and fTTl |wa(i')| ^ Cx^ which gives the following 
estimate 



t 



Jo X Jo X X 



This ends the proof of the first claim. For /x < 1/2 observe that 

POO POO 

^ / (e-''/^* - l)wxiv)dv ^ - / wx{v)dv = x^-^^^il/A - n^)/2x. 
Jo Jo 

Consequently, we get 

and this completes the proof of the second claim in the case /i < 1/2. □ 

Remark 1. If 1 < x < 2 then the absolute value of the error term is bounded by Ct^^'^, while 
for X > 2 it is bounded by C(|)^. Observe also that the density ^\t) up to a multiplicative 
constant is close to the density of the hitting distribution of 1 by the one-dimensional Brownian 
motion starting from x, when the fraction | is small with the error precisely estimated by the 
above lemma. With some additional effort one can show that the error term estimate can not 
be improved. 

Proposition 5. For every /i 7^ and c> we have 



lim 



Proof. It is enough to show the above-given convergence for strictly negative indices. The 
general statement follows from ([H]). Now we assume that /i > and consider q'x^\t). We 
define 

y,/zle~^' cosfvru) f°° e'^KJu)e~'"''y/udu 

Wc^[v] = - > -^^ —e \ t-j I fix , V 



K^^i{zi) y/2Ti Jo cos2(7r/i)i^2(M) + (7rJ^(M) + sin(7r/i)i^^(M))2 



for every f > 0. If fc^ = then the first sum is equal to zero. Using the asymptotic expansion 
we easily see that 



lim x[vj = — y — r^e 



■VZi 



The relation implies that |/u(^)| ^ C2—f= and consequently 

^/x\Ii2{xu)K^{u) - Ii,{u)K^{xu)\e~'''' ^ C2K^{u 



for every u > 0. Moreover, using the estimates of K^) and given in Preliminaries, we observe 
that the function 

f{u,v)- 



cos2(7r/i)K^(u) + {ttI^{u) + sm{iTfi)Kn{u)y 
8 



is bounded, as a function of u, by c^e (*'+2)"^3/2 |-]^^ ^^^^ ^^^ii+i/2 conse- 
quently is integrable on (0, oo). Using the dominated convergence theorem we get 

A , , cos(7r/i) p e"/s:^(M)e-''"v^c/M 



A cos TT/i r 
lim Ti7:W2\(v) = / 



^ Jq cos2(7r/i)/r2(M) + (7rJ^(M) + sin(7r/i)K^(M))^ 
which imphes 

A 



hm 



—wx{v) = Woc{v), v>0. 



Using fl33l) and fl39l) from Appendix we get that 

A, A, A, ^^,0 1 



(16) 



for some positive 6^. Next, we take advantage of the following formula of the Laplace transform 
of Wx{v) (see Lemma 3.1 in [S]) 

e-'^''wx{v)dv = , ("^ ' -{r- {j/ - 1/4)—), r > 0. 



Thus, taking the limit as x tends to oo in the above relation and applying the dominated 
convergence theorem we get 

e-..(.M.= ^-^-r + ^^ (17) 

for every r > 0. Now let c > 0. Observe that lima;/t^c,a;-5.oo ^ = cxi- Using f|T6|) and the 
dominated convergence theorem we get 

lim — — / e-'^l'^^wx{v)dv = / e-''''wooiv)dv. 



We have 

Multiplying both sides by ' tS'king limit as x/t — j- c and x — oo and using f[T7|l we get 

1 q'^-^\t)V2E /i'-l/4 , P , , , e-'^ 

lim —r-^ y— = C-- h / e WooMdv = \ — --— - > 0. 

□ 

Lemma 6. For every C > there is a constant Ci > depending on C and fi > such that 
whenever x < Ct, x > 1. 
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Proof. By Lemma [Hit is enough to consider fi > 1/2. Let < C < C. The fact that the hmit 
given in Proposition |5] exists and is strictly positive imphes that for every c G [C , C] there exist 
Ec > 0, Dc > 1 and Xc > 2 such that 



for every (x/t, x) G (c — Ec, c + Ec) x (xc, oo). The family 

{{c — Ec,C + £c)}c£[C',C] 

is an open cover of the compact set [C, C]. Consequently, there exists a finite subcover {(c^ — 
£cfc,Cfc + EcJ}k=i,...,m- Setting C* = max{xcj^ : k = 1, . . . ,m} and D = max{Dc^,k = 1, . . . ,m} 
we get 



1 1 q^^HtW27rt 

jj ^ -^^ji=rri — -pmt — ^ ^ ^^^^ 

for every C ^ x/t ^ C and x > C*, which proves the lemma for this range of x and t. Choosing 
small enough the constant C depending on /i we infer that using Lemma H] we complete the 
proof in the case x/t ^ C and x > C*. 

The estimates for x ^ C* and C ^ x/t ^ C can be deduced from the absolute continuity 
property for Bessel processes with different indices. Indeed, from ([7]) we have 



where = ^ > 0. Moreover, we have 



e 



V27rF 

and exp(— c^t) ^ exp(— c^C*/C"). This ends the proof. 



□ 



In the next lemma we show a result which provides a satisfactory estimate when t is large 
relative to x. This is done under some additional assumption on [i. 

Lemma 7. Suppose that /i — 1/2 G N. We have the following expansion 

There is a constant c > such that for t > 0, 

|E(x,t)| ^c|. 

Proof. We use the following result proved in Lemma 4.4 of [5]. Let I = fi — 1/2. Then 
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2'+i(/ + l) 



k'-'^ wx{v) dv = Ci+i{x) > 0. 



Let 



Using (|T9l) we may write 



wx[v) e 



-N./2t 



X- 



'2TXt 



■H{X,t) = X- 



/2'Kt JO 
-X^/2t poo 



wx(v e 



os;i^/+i 



wx(v e 



-K/2t 



1 / K,\i 







dv 



^2TTt 



-Ci+it 



-i-i 



i-f^){t)-X—=Ci+,t 



-i-i 



27rt 



where we applied (fT^ in the last step. Observe that 
for some constant c. Next, 



'K/2t 



1 , K 



«y «i 



which together with the estimate (see Lemma [H] in Appendix) 

\wx{v) \ = \wi^x{v)\ ^ cx^~^/2e~^'^^ 
leads to the following bound for H{X,t): 



(19) 



\H{X,t)\^ 



2t/ 



fl+2 



fl+2- 



(20) 



To complete the proof we need to find the constant Let Tq~^^ denote the hitting time of 

if we start the process from x. Due to the strong Markov property and the scaling property 
we obtain the following equality of the distributions: 



^0 "3,2^0 +^1 ) 



where Tq and are independent. It follows that 



pi-^^)^Tt^) > t) ^ P!^-^'\Tt^^ > xH) + P^-^\t[~^^ >t), t^oo. 
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Note that by the result of Getoor and Sharpe [8] we know that t^Px (Tq ^' > t) — r(^+i)2M' 
which imphes that 



x^f* - 1 



T{i2+ 1)2^' 



Fom 020]) and 



we infer 



t^'Pi->'\T[~^^ >t) = -^=Ci+i. 

/iV 27r 

This in turn shows that Ci+i{x) = '^'^^^^ly^^^^ = '^T•lJ)2^'^^ ^'^'^ completes the proof. 

□ 

Remark 2. Let / = /i — 1/2 G N and G N. Since all moments of n with respect to w\{v)dv 
exist we can write 

for some constant c, depending only on /i and k, where the constants C/+j = C;+j(x) can be 
found by similar considerations as above. 

Lemma 8. Let fi — 1/2 ^ N and let I = [/i + 1/2] . There are constants Ci, C2, C3 depending only 
on /i stic/i that 

^^^^^ ^ (1 - ^3 (y) J ^ ^i-^Ht) ^ c.^^e-' 

for t > X > 1. Note that / - /i + 1/2 > 0. 
Proof. Applying ( [T4|) we have 



„-A^/2t /"OO / 1 , 1,. . A 



The upper estimate of q^f^ (t) is obtained using almost the same arguments as in the proof of 
fl2Ul) in Lemma [71 The resulting bound is of the following form: 

kiTHt)|^cx^-^A^f^) . (21) 



t^+i V t 
12 



1 , 



Next, we deal with q[.2\t). Observing that 

-'-E/-)'^ ,,,, 

and using Lemma [16] from Appendix, where the estimate of w\^2{v) is provided, we have 
We need to effectively estimate the integral 

Using the folowing change of variables j = s we obtain ^^^-^dv = ds, which yields 

^ ^ ^ ^ 2^. (22) 

V S V 

2 



Assume that x > 2. Thus, k = 2Xv + v ^ v{v + x) and 



^ X- 



J(t,x)= / ^ ^)^+3/2 ^^^rft; = Jl(t,x) + J2(t,x), 



where 







and 

^2 ^,3; = / '---rr, -, -dv ^ x^^^^ \ — -dv. 

Applying the above change of variables and (l22l) we obtain 



J2(t,x) ^ x^^-i / -dv ^ x^^~' / ^ —ds = x^^-'^TT^ / ds. 



Observing that s+i '^^ ds < oo we arrive at 



Mt,x)^x'^-'-^, t>x. 
The first integral Ji{t,x) we trivially estimate 

Jlit,x) < CX^^-' j^^y+'/'~^^^-l/2-A^ ^ Cj2(t,x). 

This yields the following estimate 



J{t,x)^x'^-'-^, t>x. (23) 
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Next, assume that x ^ 2. Thus, k ^ v , v ^ 1, and 

/•°° 1 (\vy+^ + 

where 

Ji[t,x)= / - — -|Nn„,o 7 — — ;t dv ^ / --— — dt; 



Clearly 

, , , /-^ 1 (Aw)'+i + w^a+i) 1 
^3(t,x) = I dv ^ c- ^ cUt^x) 

if t > 1. Obviously this implies that fl25]) holds in the case x ^ 2. 
Using this estimate we finally obtain that 

- Ax^^"^ ^^^^ , t > a; > 1. (24) 

A combination of fl^T]) and f l24p ends the proof. □ 
Lemma 9. Let = 0. For t > 2a;, 

1 + logx 



X t (1 +log|)(l + logt)' 



Proof. Recalling the representation f[TU|) for gL (t) and observing that e — 1 ^ 
have (note that Wxiy) ^ 0) 



/27rt 

\-^^[x-^'^/2t+ -{-wx{v))dv 



V Jo (f^ + i) 

Hence, it is enough to estimate the integral 



J(t,x)= - — —{-wx{v))dv. 
Jo ('^ + 1) 

We start with the case x > 2. Due to Lemma fT7l (see Appendix) we have 

v< 3/2, 



> X > 2. 



We write 



J(t,x) = 




{-wx{v))dv 



= Ji{t,x) + J2(t,x) + J3(t,x). 
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Note that 
which imphes 

/•^/^ 1 K r^'"^ 1 vx 3 
^ r 1 _ r 1 vx 

J3/2 log V {k + t) 73/2 log f (f X + t) 

^ , , , , /""^ logx , [°° \ogx v"^ , 

J3{t,x){t,x) ^ / — - — —- — —rdv^ / — - — — dv. 

Jx XV^ \og^ v{K + t) XV^ log'' V [V^ + t) 

Assume that 2x < t < x^ . First, we deal with 

1 z"*/"^ 1 1 1 r 1 1 

J2(t,x) ~ — TTTT / , rdv -\ -7- / , -dv 

X ' Jz/2 f ' logf(fX + t) X^/^ Jt/^ f ^Z"' logf (f X + t) 

1 r*/^ 1 , 1 r 1 , 

-at' H — / ^TT^r^ dv 



x^/H 73/2 f^/^logw x'^l"^ J^|^ v'^l'^logv 

1 r*/^ 1 , 1 r 1 , 

-ctt> H — -7- / ^TT— dv. 



We have 



and 



X^/H J 3/2 V^/'^logV X'^/^ y^/a; f log 



dv ^ Vt/x—^ 



3/2 V^/^hgV log I 



which shows that 



1 , ^— 1 
^— dt; ^ 2^x t- J 

4/^.^3/2 logy logi 

1 1 



xv^ log 7 



Next, 



T / \/ \ I logx 1 , logx , ^ , 

J3(t,x)(t,x) ^ / — ; — 2~ / 9 , / — 0— 



X log^ V {vP- + t) a; y^, t;2 log^ V x^ log X 

Combining all the estimates we see that for 2x < t < x^ we have 

J(t,x) ^ J2(t,x) ^ \ ■ 

xVt log - 

Next, we assume that t > x^. 

V. ^ ~ /"^ logx 1 , ^ r log^ 1 . 

J3 t, xjlt, x) / 5—^-;; rdf + / ^^^^ ^dv 

^ ^^ ^ A xlog27;(^2 + i) xlog^;(^;2 + ^) 

logx /"^ 1 , logx /""^ 1 

-at; H / — T^dv. 



xt log^v X f^log^f 
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Observe that 



log^ V 



dv ^ 



1 , Vi 
-dv 



log V log t 



and 



As a consequence we obtain 



-dv 



^log^v Vtlog^t' 



Next, 



J2{t,x) 



1 logx 

Xy/t log^ t 

1 1 

X-*^/^ 73/2 ■y-'^/^ logf {VX + t) 

1 1 



J^{t,x){t,x) 



1 



(it" 



a^^'^^^i^^ ./3/2 ^'^^^ log w 



(if 



tlogx 



^ C 



Recall that 



1 log X 
x\fl log^ t' 

3 



x^ ^ t. 



Hence, in this case, it is easily seen that the integral J^itjx) dominates and 

1 logx 



J{t,x) ~ J^{t,x){t,x) 



^ t. 



x^/t log^t' 

Summarizing all the estimates obtained for J(t, x) we have that for 4 < 2x < t < 



X- 





1 1 



x-'/y2t + 



while for t > x^, 



^27rt V ' xv^logf 

Ae~^^/^* 1 + logx 

X t (l + log|)(l + logt)' 



V27rt 



X 



^/V2t+ / {e^""^^' - 1) wxiv)dv 
1 logx 



x-i/V2t + 



X 



log^ t 



/27rt 

^ log^t 

Ae"^^/^* 1 + logx 

X t (l + log|)(l + logt)' 
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This completes the proof in the case x > 2. 

Finally, assume that 1 < x ^2 ^t. By Lemma fTTl (see Appendix) we have 



1, v^2, 
V >2. 



1 



t)2(log^ V+l) ' 

Thus, using the fact that for such range of x we have 



we obtain 



^Ht) - -^rf. 



v/t yx^H Jo v^ + t A (t;2 + t)(log2t; + l) 
Obviously, we have 



2 ^2 



-dv 



v^ + t t 
and 

d'u 1 1 ds 



2 (t;2 + t)(log2t; + l) 7^7 2 (s2 + 1) (log(s) + log v^)2 + 1 ^/ilogH' 
Hence, for 1 < x ^ 2 ^ t, we have 

?x W ~ A- 2- ~ A 2 

Ae"'^^/^* 1 + logx 



X t (l + log|)(l + logt)' 
The proof is completed. □ 



4 Applications 

4.1 Survival probabilities of killed Bessel process 

In this subsection we introduce uniform estimates for the survival probabilities of Bessel process 
killed when exiting the half-line (l,oo). The theorem below is formulated for processes with 
non-negative indices, however due to ([H]) we can easily derive the corresponding result for strictly 
negative indices. 

Theorem 10. Let /i > 0. Then, for every t ^ and x > 1, we have 

Pi^^tKTi"^ <oo) "'"^ ^ 



VxAt + X - Itf" + x^^'' 
Moreover, for every t ^ and x > 1, we have 

p(0).y(0) X ^ log^ 

[1, >t) ^ ^^log(l+tV2)- 
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Proof. Using f[T^ we get 

/CO 

= Ji(t, x) + J2(i, a;). 
The integral J2{t, x) can be estimated as follows 



JtWx 

{x -ly-^i" f {x -ifi" ^\ _ X -1 / 1 ^ 



X \{t\/xy J X \{tyxy' {x-iy^ 

x-l 1 



X {ty xy + {x - ly^'' 

Observe also that for t ^ x we have 

1 1 1 



{t W xY + {x - lyf" + t^^ + x^^''' 

X — 1 X — 1 X — 1 



X ^/x-l+x VxTvt + X - 1 

The fact that for t ^ x the integral Ji{t,x) vanishes together with the above-given estimates 
of J2{t,x) end the proof in that case. If t < x, then 

X— 1 1 X— 1 1 

J2[t,x) ^ 



X x^' + {x- 1)2/^ X + x^/^ 
Substituting m = (x — l)^/(2t) we can rewrite Ji(t,x) in the following way 

x-l r _r^-n2/9, , ^2 



^ ' a;M+l/2 / a;M+l/2 /(^„i)2 



For X ^ 2 we have 

and it means that Ji(t,x) is dominated by J2(t,x) in that region. Consequently, we obtain 



Pi^\t < Ti^'^ <oo) ^ Ut, 



X — 1 1 X — 1 1 



X 



whenever t < x and x ^ 2. Moreover, for l<x<2, t<x and (x — 1) > \/t we get 

•1/2 foo 
'1/4 



1 /•1/2 /-oo 
2'^ A/4 ^ ' J, 
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Thus, the integral Ji{t,x) dominates J2{t,x) and we have 



Finally, for 1 < x < 2, t < x and (x — 1) ^ \^ we get 

Ji{t, x) ^ ^ u^^^^du = V2{x - 1) (-^ - -^J , J2{t, x) ^ X - 1. 

Thus we get 1 — 1/a/x ^ (a/x — 1)/Vi ^ V^/Vi and consequently 

1 1 \ X — 1 X — 1 



Vt Vi + x - lt^' + x^^'' 



Now we deal with the case /i = 0. We begin with the case of large time t ^ 2. We have to 
consider three cases. For s ^ t ^ x^ we have 

xs log s 

Consequently 

^ X I .log2. xlogt ^log(l + tV2)- ^^^^ 

If 2 ^ t ^ x^ and additionally x ^ 2, using the above estimate, we have 

1 > Pf (Tf > t) ^ Pf (Tf > x^) ^ -J;^ ^ 1. 

log(l + x) 

Finally, for 2 ^ t ^ x^ with x < 2 we get 2 < t < 4 and we can write 

pf (T^r >t) = < ^ + pf ^ 10) ^ A. 

To justify the last approximation observe that, using Lemma HI we get 

g(°)(s)~A, 2^t^s^l0, 

and it gives 

/iU 

Moreover, using (12^ we also get 

Pf (T^°^ ^ 10) ^ A. 

Combining all cases we obtain that 

>tj ^ ^^log(l + tV2)' 



"10 

is ^ A. 
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In the case of small times t ^ 2 and 1 < x < 2, by Lemma HI we have 
and thus 



Pf(t ^10)- I x'-^ds ^lA^^ 



Observe also that Pi°H^i°^ ^ 10) ^ A by our previous estimates in the case of large times. 
Hence 

pi°)(Tr>t) - 

Finally, for t ^ 2 and x > 2, using the Markov property, one can easily obtain that 

Pf(T(°)>t)^pf(Tr>2)^l. 
Again combining all the cases we easily obtain 

pio)(r^o) . 1 A ^ ^1 A ^og^ ^<2 

t-, [1, >t) ~ -^/\log(l + il/2)' 

This ends the proof. □ 

4.2 Poisson kernel for hyperbolic Brownian motion with drift 

Let us consider a half-space model of n-dimensional real hyperbolic space 

e'^ = {(yi,...,y„_i,y„)GM":2/„>0} 

with Riemannian metric 



ds 



2 



dyf + ... + dyl_^ + dyl 

/2 



Vn 

The hyperbolic distance d^-n{ii^z) is given by 

|2 

22/n^n 

The Laplace-Beltrami operator associated with the metric is given by 



\y — ^\ 

cosh den (y, z) = 1 H — , y,z^W^. 



^ dyf dy„ 



1=1 



For every /i > 0, we also define the following operator 
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The hyperbolic Brownian motion (HBM) with drift is a diffusion F^^) = {Yi^\t ^ 0} on W 
with a generator | A^. For /x = we obtain the standard HBM on (with ^ Aen as a 
generator). 

The structure of the process Y^^'^ starting from (y, G H" can be described in terms of 
the geometric Brownian motion and integral functional A^^^^ as follows. Let B = {Bt,t ^ 0} 
be {n — l)-dimensional Brownian motion starting from y G M""^ independent from a geometric 
Brownian motion X^~^^ starting from ?/„ > 0. Then we have 

Yl'^^iBiA^ymxt'^). (26) 

We consider D = {{yi, . . . , yn-i, Vn) G H" : yn > 1} and the first exit time of Y^^^ from D 

td = inf{t ^ : f/^^ ^ I^} = inf{t ^ : X^^"''^ ^ (1, oo)} = r, 

where r = is the first exit time from the set (l,oo) of a geometric Brownian motion X^"^) 
defined in Preliminaries. We denote by P^^\y, z), y E D and z G dD, the Poisson kernel of D, 
i.e. the density of the distribution of Yf^, with Yq^^ = > 1. 

Theorem 11. For every fi > we have 

P^'^ {y.z)^ { J\ ''''' , (27) 

where y = {y, y^), yn > 1 and z = {z,l), z e IR""^ 
Proof. Let us denote by 

, , exp(—\w\'^/2t) ^„ -I 

(27rt)("-i)/2 ' 

the Brownian motion transition density in H.""^, n = 2,3,.... Using the fact that B and 
A^^\t) are independent we obtain 



poo 

P('^\y,z)= gt{~z-y)q^-f{t)dt, 2/„>l,^G 

^0 



Using the estimates given in f lT^ we obtain 

P^^\y,z) ^ X / e"(l^-Sl'+^')/2* ^ 



tM-l/2 + ^/^-l/2 t(n+2)/2 

J'" " - dt 



f{n+2)/2 ' yn J - ^(n+l)/2+/i 



^"/2-lg-«^^ + ^l/2-MW2 / yn/^~'"^+>' du 



= \ ,,/^-"/2-l/2 
|_g _ y|2 _|_ _)^2 

where p = , A = ?/„ — L Using fl2Sl) (see Appendix) we can see that 



M-l/2 9 
^ (|5-y|2 + A2)-/2' |5-^|2 + A2^ 
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and 



p(.)(,,).A (^^")^"^^^ <i 

Combining both estimates and using the formula for the hyperbohc distance we obtain 



□ 



Remark 3. The operator is strongly elliptic operator on every bounded (in hyperbolic metric) 
subset of H". Consequently, the hyperbolic Poisson kernels of such set are comparable with 
Euclidean ones. However, considered set D is unbounded in and the general comparison 
results can not be applied. Besides, the function P^^\y, 2) for /i 7^ 1/2 is no longer comparable 
with Euclidean Poisson kernel of upper half-space and the difference in behavior of those two 
functions is described by the factor 

M-1/2 



cosho?H"(?/, z) 



5 APPENDIX 

5.1 Uniform estimates for some class of integrals 

Lemma 12. For u^O,0^a<b and d > we have 



u^e-'^du ~ ( ^ e-- \J " ■ (2J 



b + y d{b-a) + l'' 



where c = c(z/). 



Proof. Let a, 6, (i) = Jl^u'^e '^du. Since F{h',a,b,d) = d ^F{h',ad,bd,l) it is enough to 
prove the lemma for d = 1. Assume that 6^1. Then 

b rb—a pb—a 

u'^e-'^du = e-" / (a + uYe-'^du ^ e'" / {a" + u'')e-''du ^ e-^ia" + 1)((6 - a) A 1), 
Jo Jo 

which is an equivalent form of f l28p in the case b'^d = l. Ub<l then 

u''e-''du^b''{b^a), 

which is exactly f l28|) in the case b < d = 1. Note that in all comparisons above the constant c 
is dependent only on u. 

□ 
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Lemma 13. Let ^ a ^ 1. Then for every v > we have 



e ^^udu 



lo log^M+l {v + l/ay(log\v + l/a) + l)' 
If additionally av ^ 1 then 

1 g-™^(^^ 1 — a 



„ 1-logu (w + 1)3/2(1 + log(i; + l))' 
Proof. Let 

, , ;° e-^'^'udu 
J[a, V 







log^ u + 1 



First, assume that av < 2 then 



2 



J(a,v) ^ / ; — 5 ^ „ 

io log^M + l 2(log2a + l) 

If additionally a > 1/2 then 

J(a, v) ^ 1. 

If a ^ 1/2 then 

J(a, V) ^ e I - — 2 ^ c — 2 



a/2 log n + 1 log a + 1 

which ends the proof fl2I?l) in the case av < 2. 

We assume now that av ^ 2. Observe that for every g, r G M, from the fact that 
(gr/Vr^ + 1 — + 1)^ ^ we get 



(g-r)2 + l^ ^ 



r2 + 1 
This implies that 

1 log^ s + 1 



(log — log s)2 + 1 log^ V 

Consequently, we obtain 



av 



J{a,v) = ^l ^ 11 ^ 2/2 / e 's(log^s + l)rfs. 

v^ Jo (logt; - logs)2 + 1 v^log^vJo 



On the other hand we get 



J[a,V) = — I — ; r- ^ ; / 6 SdS, 

v^ Jo {\ogv-\ogsY + l^ v-'ilog^v + l) 



which ends the proof fl29|) . 
Let 



j( \ i' ^ 
J{a,v) = / du 



a 



1 — log U 

ds. 



w3/2 J^^ 1 — log S + log V 

23 



If f < 2 then 



J(a, v) 



1 — loer u 



du ^ 1 — a, 



(31) 



which is fl30|) in this case. 

Next, we assume that 2 ^ v ^ 1/a. Using the fact that for every r > g > we have 
r — g + l^r/(g + l)we get 



1 - log s + log W ^ \Jq J I J 1 - log s + log V 



ds 



log V 



I e-'s^'^ds+-^ I e-' s^'^ (1 + log s)dsf^ 
Jo logwJi 



Moreover, for 2 <v <l/avfe get 



^-.^1/2 



1 — log s + log V 



■ ds ^ / ^ ^ as ^ 



We have just proved that 



J(a, v) 



■y3/2 log?; 

Combining (EI]) and (152]) we complete the proof of (|5Uj) . 



1 1 — log s + log V log w 



2 < < 1/a. 



C2 

log V 



(32) 
□ 



5.2 Estimates of wi^x and i(;2,a 

Lemma 14. There exist constants c = c{fi) > and 6^ > such that 

\wi xiv)\ ^cx'^-^^'^e-''^^, v>0. 



(33) 



Proof. Recall that the set of zeros of the function K^{z) is denoted hy Z = {zi, ...,Zk^}. For 
every Zi E Z we have ^Zi < and the set Z is finite. Consequently, there exists constant ci > 
such that for every i = 1, . . . , we have 



^ ci(x - 1). 



Kf,{xzi) 




K^{xzi) 


- K,{z,) 


K^,-i{zi) 









Moreover, using ([S]), the constant Ci can be chosen to ensure that for every x ^ 2 we have 



e^^^K^{xzi) 



Ka-l(Zi 



X 



Hence there are C2 > and = —maXi{^Zi} > such that 



ZiC^'^'KJxZi 



^ 



X 



3/2 



and it gives 



□ 
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Let us define for n > and x > 1: 

S^{x, u) = (xu) KM - (^") • 

The function appears in the formula for W2,a and consequently, the uniform estimates of the 
function S*^, which are given in the next Lemma, are crucial to get the estimates of the function 
W2,x given in Lemma [16] for /i > and in Lemma fT71 in the case yU = 0. 

Lemma 15. For fi ^ we have 

— , ^ ^ bn[x,u) ^ A—— — -. (34) 
X KM ' ^ ^ KM) 

There are constants Ci and such that for ^ 0, 1 < x < 2 and u > we thus obtain 

CiAe"^" ^ ^ CaAe^". (35) 

There is a constant Ci such that for fi > 0, x > 2 and u > we have 

cilf, (xu) K^{u) ^ S^{x, u) ^ {xu) Kf,{u) . (36) 
There is a constant ci such that for fi = 0, x > 2 and xu > 1 we have 

ci/o (xu) A'o(u) ^ So{x, u) ^ Jo {xu) Kq{u) . (37) 
For fi = 0, X > 2 and xu < 1 we have 

So{x, u) ^ logx. (38) 



Proof. Write for m > 



Then by ([T]) we have 



ip{u) 



JM. 

KM 



1 1 

u KM ' 



Writing 

Sf,{x,u) = [4j{xu) - ij{u)] Kf,{xu) K^{u) 
we obtain from the Lagrange theorem 

xu — u K,,(xu) K,,(u) 



Sf,{x,u) 



6xu K^{9xu) K^{9xu) 



The quantity 6 here has the property 1 ^ 6x ^ x. This and the monotonicity of the function 
Kfj,, give the estimate (134|) . The estimate (135|) is a direct consequence of the hmiting behaviour 



of the function K^. 



To prove fl5Bl) and (1571) note that the function g{x, u) = = k^^u^/^(xu) ^ function of 
X is decreasing. Hence for x > 2, 

g{x, u) ^ g{2, u) < l,u > 0. 
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If /U > 0, then the hmits at and oo of g{2, u) are strictly less then 1. By continuity 

sup (7(2, u) = a < 1. 



u>0 

If /i = by the same argument for x > 2, 



g{x,u) ^ sup g{2,v) <l,u^ 1/2. 



If > 1, M < 1/2 and x > 2 then 



/oHA^o(xtx) /o(l/2)Jio(l) 
Kq[u)Iq[xu) Ao(1/2)/o(1) 



These estimates imply that 

g{x,u) ^ sup g{2,v) = a < l,u > 1/x. 
Hence, in both cases (/i = or /i > 0) we have for x > 2 and xm > 1 



Si,{x,u] 



1 - 



ip{xu) 



4}{xu)K^{xu) K^(u) ^ (1 - a)I^{xu) K^{u) 



which ends the proof of (136|) and (137|) . 

It remains to consider /i = 0,x > 2 and xm < 1. We apply the asymptotics of Kq at 0. 
Namely, by (|4]) we can write 

K^{z) = -\og^-Io{z) + A{z), 

where A{z) — )■ c > 0, 2; — )■ 0. This yields 

5'o(x, u) = Iq{ux)Iq{u) logx + Iq{ux)A{u) — Iq{u)A{ux) ~ log 

Lemma 16. For /i > and x > 1 we have 



X. 



□ 



^2/^-1 

W2Av) ^ (- cos(7r/i)) ^^^^^^^3^^^^^^^^^^^^ , ^ > 0. (39) 



Proof. Let us denote 



h{x,u,v) = '^^ '^^^ e"^"e"''''u 

cos^{7i fi)K'^{u) + {ttI^{u) + sin(7r/i)A'^(M))2 



Since W2,x{v) = — cos{nfi)'^ h{x, u, v)du, it is enough to estimate h{x, u, v)du. which is 
done below for two cases. 
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A) Case 1 < x < 2. 
Suppose that < m < 1. By Lemma fT5| 5^ {x,u) ~ A, thus 

^ AM2^+^e-(''+^)", < M < 1. 
For M > 1 we have, by Lemma [T^ {x,u) ^ cAe", which yields the following upper bound. 

h{x,u,v) ^ c e-^"e— u 

^ cAMe-(''+^)". 
Applying Lemma [12] to the above estimates we arrive at 

h{x, M, v)du 



and 

r-oo \ 

u, v)du ^ c- rTTe"(''+^\ 

1 + 1)^ 

Combining both integrals we obtain 

h{x, u, v)du ^ - ^ 







which proves the lemma in the case 1 < x < 2. 

B) Case x > 2. 
By Lemma fT5| S'^(x,m) ~ I^{xu)K^ (u), which implies 

iq{u) + i^{u) 

Next, using the asymptotics of the Bessel functions, we arrive at 

^-l/2^/^+l/2g-t;«^ I/X < U < 1, 

To estimate H{x,v) = h{x,u,v)du we split the integral into three parts: 

rl/x p1 poo 

H{x,v) = I h{x,u,v)du+ I h{x,u,v)du+ / h{x,u,v) 

Jo Jl/x Jl 

= Ji{x,v) + J2{x,v) + J3{x,v). 

Applying Lemma [T^ with a = 0,b = 1/x and d = v + X, the first integral can be estimated in 
the following way: 

Ji{x,v) ^ x^" f ^ M^'^+^e-^'^+^^^rfM 
Jo 

x^ 1/x 
{v + A)2m+i 1 + {v + \)/x 
x^ 
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Next, we deal with the second integral. Again, by Lemma ITU with a = l/x,b = 1 and d 
we obtain 

J2{x,v) ^ [ 

Jl/x 

1 + i; — f /x 

fl/2 ^ 





1 




V 


1 + 


1 


V 


^+ 


1 




V 


1 + 


1 




V 


1 + 


V 
X 


1 + 


V 



X + \v 

1 



+ 1 

The third integral can be estimated for f > as follows. 

POO 



(^ + 2) 



e 



-{v+2) 



It is clear that 

Mx,v) ^ c 



+ l)A«+3/2(^ + 2,)M+l/2- 

Next, 



+ l)/^+3/2(^ + _ + + X^ _{V + 1)^+3/2 

Xf" ^ ^ Xf" {V + x)2m+2 ~ [v + xY+^l'^ 



and 



(,,+ l)A>+3/2(^^^)M+l/2 ^ _ {V + + a:)^+V2 ^^^^ / 1 + £ X ^ 

^ xi" ^ \l + v) ^ v + l 

= {l+v/xf^+^e-"/''. 

The observation |^|||p^p+3/2 + (1 + v/x^^'^'^e''"^^ ~ 1 completes the proof. 



Lemma 17. For = and x > 1 we have 



, , 1 log(x + 1) 

-w^fw) — — ■ — > 

^ x(w + l)3/2(^; + x)V2log(y + 2)(log(x + l) + log(t; + 2))' 

Proof. Let us denote 

1 log(x + 1) 



f{x,v) 



x{v + l)3/2(t, + x)V2 iog(t; + 2)(log(x + 1) + log(t; + 2)) ■ 

We write 



Kl{u) + Ti^mu) \J, Jy, J, ) Kl{u) + 7i^P,{u 

= Ji{x,v) + J2{x,v) + J3{x,v). 
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The estimates of Jslx, v) are exactly the same as the corresponding estimates proved in Lemma 
fT6l Hence, for x ^ 2, 

1 e^(''+2) 

and for 1 < X < 2, 

Mx,v) ^ c^^^e-^^^'l (41) 
To estimate J2{x, v) for x > 2 observe that for xu ^ 1 and n < 1, So{x, u) ^ Kq{u) ^ 



xu 

(1 — logw). This follows from the asymptotic expansions for Iq and Ko (see (H^ and (l3])), 

/xu 

and Lemma [151 Thus 

xl/2(^;+ 1)3/2(1 +log(y + l))' ^ ' ^ ^ 

where the last step is a consequence of fl30l) with a = 1/x. Next, 



11 

^ c— 7:7- ^,„,„,^ ^ — 7 —7, u > 0, (43) 

(t; + 1)3/2(1 + log(v + 1))' ^ ^ 

where, again, the last step is a consequence of (l30l) with a = 0. 

For X < 2 and 1/x ^ m < 1, by Lemma [151 we have S'o(x, m) ~ A hence 

1 g-™^ /.I g-™^ 



J2ix,v) ^ X — 97-2/ ^ du^ X / ; — —du. (44) 

A/, A2(n) + 7r2/2(u) ii/, l + log^M 

Finally, for < u <l/x we have, by Lemma [T51 5'o(x,u) ~ logx, a; > and consequently 
using (129|) . with a = l/x, we obtain 



Ji(x,w) f« logx / , „ ^„ , , eiw ^ log X / ; — 7y—du. (45) 







logx 



(t; + x)2(log''(t; + X) + 1) 



(46) 



We are now ready to estimate the function —wx{v). At first we consider 1 < x < 2. Taking 
into account that logx ~ A, using (l44l) and (l45l) w arrive at 

Ji(x, f ) + J2(x, f ) ^ X 2 — du (47) 

Jo 1 + log u 

' (48) 



(t; + l)2(log'(t; + l) + l) 
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Combining this with (HT]) we have for 1 < x < 2, 



-Wx{v) = j{Ji{x,v) + J2{X,V) + J3{X,V)) ^ j{Ji{x,v) + J2{X,V)) 

1 

f{x,v). 



(t; + l)2(log'(t; + l) + l) 

Taking into account (HO]) . and (jlSD, we infer that for < < 2 and x > 2, Ji{x,v) + 
J2{x,v) ^ cJ3(x, f), which yields 

-w^A(f) = jiM^^ + <^2(a:, t') + Jsix, v)) ^ ^ Jsla;, v) ^ ^ /(x, v). 
For 2 ^ < X, by (gnD, (gOD and Ji(x, v) + J3(x, i;) ^ cJ2{x, v). Hence, 

-wa(w) = \{Ji{x,v) + J2(x,t;) + J3(a;,t;)) ~ \ j2ix,v) ^ ^ ^ f{x,v). 

Finally for t> ^ x > 2, by (l43i) . (146|) and (HOl) we have J2{x,v) + J3(x,f) ^ cJi(x,f) for some 
constant c > 1 and consequently 

1 1 log X 

= -(Ji(x,t;) + J2(a;,^^) + Jsix.v)) ^ -Ji(x,t;) ^ — - — — ^ /(x,w). 

A A xw log f 

The proof is completed. □ 
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